» NASA CE- 1323 99 | 



Mechanical Engineering Department 
South Dakota State University 
Brookings, S. Dak. 57006 


PRIOR-TO-FAILURE EXTENSION 
OF FLAWS UNDER MONOTONIC 
AND PULSATING LOADINGS 


Inelastic Fatigue 


Michael P. Wnuk 


(if AS A^rp^. * 

S5 ^ 3 fla «s To 'Failure — — 

stat e aI r Pro «e s C ff : Urif s °g*°Mc UB 
ate V*iv.)&%* e Port (sourfc f ^Igue 
I BC W.00 Daicot a 

I ' CS CI 20K 



2927 j S 



National Aeronautics 
and Space Administration 

Grant NGR 42-003-006 

Annual Progress Report 

May 1971 



Annual Progress Report 


PRIOR-TO- FAILURE EXTENSION OF FLAWS 
UNDER MONOTONIC AND PULSATING LOADINGS 
(Inelastic Fatigue) 


May 19 71 


This work was supported by 
National Aeronautics and Space Administration 
Grant NGR 42-003-006 


by 

Michael P. Wnuk 


Department of Mechanical Engineering 
South Dakota State University 
Brookings, South Dakota 
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by 

Michael P. Wnuk 
(Abstract) 


An equation governing the prior to failure crack propagation is proposed. 

For a rate-sensitive solid containing two-dimensional crack and subject to the - 
tensile mode of fracture the differential equations are integrated numerically 
for the loads increasing monotonically in time. The resulting integral curves 
a = a(£) and & = £(t)» i.e. load vs. crack length and length vs. time, indicate 
that the growth of cracks in the subcritical range is strongly rate dependent. 

The fatigue growth, viewed as a sequence of slow growth periods, is simu- 
lated on EAI 380 analogue computer. The fourth power law proposed by Paris is 
confirmed only within certain range of high- cycle fatigue propagation and for a 
rate-insensitive solid. Otherwise, that is for a more pronounced rate dependency 
induced by viscosity of a solid and/or in the proximity of the final instability 
point the growth is markedly enhanced. For sufficiently small ratios of the 
applied stress intensity range AK to the toughness K c , the suggested fatigue 
growth law consists of two terms, i.e. 


di 

dn 


■ if {«£>* * «-'#»■] ■ 


= ttK 2 /8Y 2 
c 


First term is the familiar Paris expression while the second one accounts for the 
rate-dependent contribution; f denotes frequency and Y is the yield strength. 


H 
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Prior- to- Failure Extension of Flaws Under 
Monotonic and Pulsating Loadings 

Part I Basic equations. Monotonic loads . 

The catastrophic fracture is often preceded by a quasi-static extension of 
an initial defect which is too small to be detected. The process of slow propa- 
gation occurs at loads below the critical level, and it connot be described in the 
framework of Grif fith-Irwin theory of fracture. In fact the differential equation 
which governs the subcritical growth in a quasi-brittle solid, [6] 

M(o,£,da/d£) + G(a,£) = G c (1.1) 

implies possible extension of a pre-existing flaw at incredibly small initial flaw- 
sizes (or, equivalently, at very low stress levels). It is only at the end of the 
slow propagation stage, when the "slow growth operator > vanishes^ that the 

Irwin criticality condition is satisfied, i.e. 

G(a ,£) = G c (1.2) 

or for the linear range of fracture mechanics 

K(a,£) = K c (1.2a) 

where G denotes the energy release rate, G is the specific fracture energy, K(a,£) 
denotes the stress intensity factor and K c is the fracture toughness. During an 
infinitesimal growth the applied stress a and the corresponding crack length £ 
undergo the change 

(a ,£) (a + da, £ + d£) (1.3) 

This is associated with the plastic energy dissipation absorbed within the end-section 
of a progressing crack (l < x < a) i.e. 


SU = (SU) £ + (6U) ct 


(1.4) 
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in where 


(6U) - 4Y J u(x,a,£) dx 

& = const d& ^ 


<«Wa = const “ 4Y u(tl P } + 2Y |l } u(x ’°’*> dx 


(1.4a) 


Here u(x,a,£) denotes the displacement of a Dugdale crack evaluated within the 
process zone. The first expression in eqs. (4) is identified as twice the slow 
growth operator, while the second one is twice the energy release rate. Require- 
ment that the energy balance is satisfied at every instant of the slow propagation 
stage leads to the governing equation (1); spell out as "Michael plus George equals 
critical (or crazy) George". The resulting differential equation is valid for 
the sub critical range of loads 


a o 1 a 1 a c * or 


(1.5) 


K 0 1 K(a,£) < K c 

in which Co (° r K Q ) is the propagation threshold, while a c (or K c ) is Irwin's critical 
threshold identified here with the transition to fast propagation. For many ductile 
solids with well-defined flat "yield shelf" the ratio of the initiation to rapid 
propagation threshold, K 0 /K c , can be estimated as 


Ko (X)^ ~ 1 (plane stress) CU6) 

K c E 30 

where Y is the yield strength and E is the Young modulus. For strain-hardening 
materials and under high triaxial constraints the above ratio may approach unity, 
which implies a neglegible amount of slow growth. 
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Examples of application . 

For a two-dimensional Dugdale crack equation (1) reduces to 

JL £ 2 [8sec 2 3 “ tanf3] + c[£tan$ 4- logcosft] = 1 (1.7) 

2 dc 

while for a penny-sh aped crack with an associated Dugdale-type plastic zone it reads 


1 c 2 x 3 Ax l - / l-x 2 _ L 

3 (l-x 2 ) 3 / 2 dt /w? 

Here, the applied load a and the crack length Z are normalized as follows 


( 1 . 8 ) 


( B in(l. 7) 

i 

i, YX in(l. 8) 


ttK| 

- C(£*) = C(— 
8Y 2 


(1.9) 


For small scale yielding range the above expressions can be expanded into a 
McLaurin series at 3 0 (or X -*■ 0) , yielding 

— = » plane crack 

dc 2 p 

( 1 . 10 ) 


Neither of these two equations could have been integrated in a closed form, but 
as we proceed to show, the numerical (IBM- 360) or analogue computer (EAI 380) 
integration allows one to derive certain simplified rules to be used for prediction 
of the sub critical growth under 

(a) monotonic loadings , 

(b) pulsating loadings. 

Within the linear domain of fracture mechanics, when yielding is confined to a 
narrcw zone small compared with crack length (q^Y ^ important generalization 
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of eqs . (10) is possible. In such a representation a viscous behavior, in addition 
to the small scale yielding already represented by equations (10) is accounted for. In 
this range it is sufficient [G] to replace the time-independent K-factor, K(cr,i), 
by the "effective" time- depen dent K-factor, namely 


K eff = K(a,£) ifr (A/I) 


( 1 . 11 ) 


where ^ denotes the normalized creep compliance I(t)/l(o), £ is the growth rate 
and A denotes the intrinsic opening distance (a material constant). With (11) the 
equation governing the slow propagation within the subcritical range becomes 


<! C6 + C/6) 


plane crack 


£A + C/A) — = — , penny-shaped crack 


( 1 . 12 ) 


dc 


C A 


where C characterizes the rate-sensitivity of the solid 


c = (A/£*) , or 

C - A [ Vi ) t _ 0 (8Y 2 /irK 2 ) 


(1.13) 


The constant A has to be supplied by the experiment. 

Two extreme cases result from equations (12) straight forward: 

(a) rate insensitive material, in which all the dissipation can be ascribed to 
plastic tine- in dependent deformation. Then C = 0, and eqs. (12) reduces to (10). 

(b) highly rate sensitive solid, say a linear visco-elastic matrix containing a 
crack with a negleglble amount of plasticity present around crack tips. Then eqs. 


(12) degenerate into 
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or 


(e/g) 



(C/A) 


dA = 2 

dC CA Z 


= CC6 2 /2 
|| * CCA 2 /2 


(1.14) 


(1.15) 


It is readily observed that the last two expressions have one common form 


dc 

dt 


C 


r 


K 2 (°, £ ) 


(1.16) 


This can be integrated for a prescribed load history 6 = 3(t), yielding the length 
of the extending crack as a function of time, C = S(t) . Examples of such integration 
are given in the second part of this report, see Fig. 11. 

Numerical Examples 

Let us illustrate now how the above equations work. Consider a central crack 
of half length l roughly ten times greater than the size of the associated plastic 
zone £*. The initial crack length is thus 


£ = (10)£* = (10) (-- f) ' (1 - 17) 

oY 

and suppose that both the fracture toughness K c and yield strength Y are kncwn. We 
want to evaluate the load at which this crack becomes unstable, say a c . As the 
first approximation let us apply Irwin's criterion for failure 

K(a,£) = K c ( 1 . 18 ) 

For a central crack contained in a large plate the above equation reads 


a 2 ( tt£) = K 2 
c 


(1.19) 
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or, in a dimensionless form. 


S 2 C = 2 


( 1 . 20 ) 


With £ defined by eq. (17), we obtain 




0.4472 Y 


( 1 . 21 ) 


The Grif fith-Irwin criterion, therefore, predicts no growth for all loads £ less 
than 0.4472 and the rapid propagation occuring at 0 C = 0.4472 and 4 C = 10. No 
slow growth can be accounted for. 

Assume new that the solid is quasi-brittle and the propagation threshold 
K 0 is about a quarter of the transition threshold K c , say K 0 = 0.224 Kc . 
Since 


Bo Kq 



we have the starting value of the applied load 


( 1 . 22 ) 


B c = B c (K o /K c ) =0.10 (1.23) 

This together with C 0 = 10 provides the initial condition for equations (10) and (12). 

Numerical integration with the rate sensitivity Cassumed to be zero yields the 
following values of the crack length and the load at instability, say and BfJ 


C f - 11.4446 


(1.24) 


$ f = 0.41793 

Thus we conclude that the catastrophic fracture is preceded by a slew extension of 
magnitude 


M = (s f - C o ) £ * 

M - (1.44)(iyK 2 /8Y 2 ) 

The final instability occurs at the load 


(1.25) 


cf f = (0 . 418)Y 

This is 6.5% less than the load predicted by the Irwin criterion, a c = (0.447)Y. 
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For larger initial crack lengths and non- zero sensitivity C. the discrepancy 
becomes much more pronounced. Some data substantiating this point are gathered in 
Table I, while the corresponding graphs are shewn in Figs, la and lb. 


Table 1 Instability preceded by the slow growth vs. Irwin's instability (Numbers are 


generated by eq. (12) and IBM 360) . 


Initial Length 
Irwin 1 « Ins tabi li tv 

C 0 “ 10 
0 = 0.4472 

o 

r-t 

II 

o 

Kf 

C 0 = 100 
B c = 0.1414, 

C c » 100 

C Q = 1000 
$ c * 0.0447, 

■ CJ 

Cc = 1000 

Sensitivity to 

— d 




Bf 

Cf ” 

Pf 

Loading Rate 
Ratio 

C/& 

Cf 

3f 

Cf 

(£f’ denotes 

0 

11.4446 

.41793 

102.739 

.13952 

1004.06 

.04463 

the set of final 

.1 

11.4856 

.41721 

102.786 

.13950 

.1004.08 

.04463 

crack length and 
applied load at 

1 

11.8188 

.41147 

102.994 

.13935 

1004.18 

.04463 

which transition 
to rapid prop a- 

10 

14.1074 

.37665 

104.979 

.13802 

1005.26 

.04460 

gation takes 
place. 



100 

24.0228 

.28854 

117.827 

.13028 

1014.48 

.04440 


The amount of slow growth which precedes the catastrophic fracture is distinctly 
dependent on the initial crack length and ductility of the solid. We have 


tiK? 

A i = f(C 0 ,C) 


(1.26) 
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1/2 , , n (1+N)/2 

or, with K c replaced by Y(it 6) (ef/£y) 


tt2 N+^ 

Aj. - f(e 0 ,C) J- s(e f /e Y ) , n*i 


(1.27) 


The analytic form of the function f(C 0 £)is not known, but for a given and C 
f(C Q >C) can be read out from Fig. 3. The other essential material parameters are 
the structural size 6 (<§:£*) and the ratio of strain at failure to the yield 

strain ty (=Y/E). 

It is obvious that the deviation from Irwin’s theory becomes more pronounced 
for lower toughness threshold levels, for larger initial crack lengths, for enhanced 
ductility and for materials which are rate sensitive. All the factors mentioned 
contribute to the inelastic behavior of a solid. 


Fart II Fatigue in Rate Dependent Solids . 

Fatigue crack propagation may be viewed as a sequence of extensions (or steps) 
of "slow growth" type, each of which occurs while the stress increases during the 
loading cycle. Therefore the amount of growth produced during one cycle can be 
computed from eqs. (1.12) by simply Integrating both sides of the equation over 
the load range 3 min f '0 f 6 m ax and regarding the current crack length ^ roughly 
constant within the single cycle. This latter assumption may not be true for the 
final stages of fatigue life, where one observes a substantial acceleration of the 
growth pace, but it certainly is all right for the major portion of the high cycle 
fatigue life. From eq. (12) we have 


(dO 


per cycle 



Kmax 

L 


B 3 dg 

2-EB 2 



jUax g2dp 

l 2-ce 2 " 
Vn 


( 2 . 1 ) 


hence 
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2 

dz _ 2 , min , _1 ,„2 _ R 2 s + 

dn 3 18 2-Z& 2 a 3 (Bm * n Bmax> 

*= w max 


( 2 . 2 ) 


+ 


~ tt — f 3min” 

\$> l 


3max 


+ log 


(2+6 max ^)(2-& ml _ n ^2z) 

^“^max ^(^min 1 ^ 


If the propagation threshold 3 0 exceeds the minimum stress within a cycle 
the lower limit of integration in eq. (2.1) should be replaced by $ 0 . This 
may considerably alter the propagation rate, see Fig, 2. 


For the initial stage of high cycle fatigue progressing at stress intensity 
remote from the criticality point, the product £|3 2 in the denominator of the 
integrand of (2.1) can be neglected vs. 2, and we have 


(dO 


per 




(2.3) 


hence 


i ■ 12 ^ ^ t*L** 3 J 

where (|l) denotes the average rate within a cycle. 


(2.4) 
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The range A3 = 3 max - 3 min does not have to coincide with the integration limits 
3 0 3 A3, if we allow for the threshold level 3 = 3o to li e above the minimum 
stress in the loading cycle 5 see Fig. 2. 

Fig, 2 The location of the threshold level 3 0 predetermines the rate of growth. 




In particular, when one restricts the attention to a zero- to- maximum stress cycle 
(3 m i n = 0), and assumes zero threshold 3 Q = 0, eq. (2.4) simplifies as follows 


dc = 1_ 
dn 12 


C 2 (A6) lt + C(A6) 3 

6<e> 


(2.5) 


The first term here can be readily identified with the Paris expression for fatigue 
crack grcwth rate, while the second term represents an additional contribution due 
to the rate sensitivity of a solid. The first term accounts for crack extension in 

onltf 

a quasi-brittle solid with no time-effects, and thus it involves the range of K-factor J # 
The second term decreases when the rate of loading is increased; it does depend on 
the frequency. The average rate of loading within a single cycle can be related to 
the frequency as follows 


< 0 > = 2fA$ 


( 2 . 6 ) 
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With this and with (A8) 2 C replaced by 2(AK/K C ) 2 eq. (2.6) becomes 

d£ _ Ji* ,AKv 4 . CA* f -l/AK\ 2 (2.7) 

dn ~ 3 S( ' 12 

c ^ 

The first term alone gives the well known Paris law*, while the second one tends to 
increase the pace of growth, at least within the frequency range in which the mater- 
ial is rate dependent (i.e. for f 1 comparable to the characteristic relaxation 
time of a solid). 

Computer Simulation of Fatigue Growth. 

As ther is no closed form solution to the equation governing crack growth in the 
sub critical range 

~ - ( 2 . 8 ) 

dC 5B 2 ( 58 + C/S) 

We employ the analogue computer technique. The program (see Fig. 4c) has been 
arranged in such a way that the integration which starts at the initial point (t Q , 8 0 ) 
proceeds up to the point of maximum load in the cycle, then interrupts and reverses 
to the new '^Initial" position 

C 0 + AC , S 0 

where Ac denotes the amount of growth within a single cycle. The binary counter 
recorded the total number of cycles before the critical point was reached. This point 
is distinguished by zero slope dlB/dC, and it is clearly visible on all photographs 

* If the range Ak = AK/K c is not very small, the Paris law ought to be replaced by 


d i 
dn 



(log 


1 

T-(Ak) T 


(Ak) 2 } 


(2.7a) 
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which depict the final stage of fatigue life. Figures 5 and 6 show an increase of 
the propagation rate due the time-sensitivity 


C = 



t=0 


A/£* 


- ttK 2 /8Y 2 


(2.9) 


of a viscoelastic-plastic solid. In the limit of C-^0 and within the high— cycle 
range one recovers the fourth power law valid for rate-independent quasi— brittle 
solid. 

Fig, 7 shows three runs at 


0.1 

< 

B 

< 

0.14 , 

AB 

= 0.04 

0.1 

< 

3 

< 

0.17 , 

AB 

= 0.07 

0.1 

< 

B 

< 

0.20 , 

AB 

- 0.10 


thus demonstrating shift in the location of the critical point and the change of 
the rate of growth due to varying stress range. Total numbers of cycles to failure 
are gathered in Table II. Each run consists of two parts: 

1 £ C <10 

10 £ C <100 

and then the total number of counts (cycles) is obtained by summing up N(l,10) and 
N(10,100). The results are repe a tab le within the accuracy of about 8%. 

A somewhat different test is pictured in Fig. 8, where four runs are shown, all 
at the same range A3 = 0.01 but at different levels of the mean stress (see also 
Table III). The dependence of N T0T on the stress range, the mean stress and rate 


sensitivity is summarized in graphs shown in Fig 9 and 10. 

Finally, Miner's cumulative damage law is tested in a series of runs performed 
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at varying stress levels. Miner’s prediction 


( 2 . 10 ) 



- number of loading cycles imposed 
at mean stress level TT^ 

N^£ " number of cycles which would lead 
to failure if the level was 
maintained throughout the test. 


turns out to be fairly well satisfied (see Table IV) . The agreement deteriorates 
a little for larger values of the sensitivity parameter and for greater discrepancies 
between two stress levels used in the program. This trend appears to conform with 
the deviations from eq. (2.10) reported by investigators who tested Miner’s law 
in the experimental way. 

The next test pertains to the limit case of highly rate-sensitive solid in 
which viscous dissipation is dominant over the plastic work. Then equation (2.8) 
simplifies to 


dC = 
dt 



K 2 (pj) 


or 


d£ _ r 
dt " C£ * 


K1 


K 2 (a,*) 


which in turn can be written as 


» _ £* r m 2 (t)x 
X £ q l-m z (t)x 


ii c 

£ q C x-m^(t) 


( 2 . 11 ) 


( 2 . 12 ) 


Equation (a) above governs the growth of a central crack propagating through a 
uniformly loaded plate; then 


(a) 
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while the equation (b) describes the growth induced by the pair of point forces P 
applied at the center of the crack surface; then 

K(o,£) = P/OrA) 3 * (b) 

Both a, (or P) and Jl vary in time. Function x(t) = £(t)/& 0 is subject to determination ? 
while the load m(t) = a(t)/a c (or m = P(t)/P c ) is given. Both equations (2.12) are 
then programmed for an analogue computer, see Figs. 4a and 4b, and the resulting 
integral curves x = x(t) for 3 various loading regimes are shown in Fig. 11 • It 
appears that the pulsating loads (sinusoidal and trapezoidal) produce more rapid 
growth in the initial stage of propagation than the constant load maintained at the 
level coinciding with the mean stress within the cycle. Yet another test, not shown 
here, proved that the randomly pulsating tensile stress produced most severe increase 
of the flaw size. 
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Appendix . Equation of Motion in Ihe Subcritical Range 
The rate of work done in separating two surfaces 


W = 


T i u i dS 


AS 


T ^ = traction vector 
u^ = displacement vector 


(A. 1) 


can be computed from a Dugdale-type model of a plastic zone at the crack tip 
embedded in a linear viscoelastic solid. If the process zone extends between 
x = £ and x = a, Y denotes the yield stress and* 

a = stress applied 

u = u(x,a(£) , £(t)) 2£ = crack length (A. 2) 

t = time 


is the displacement of the crack face, then (A.l) takes the -form 


W = 4 J Y • u(x,a,£)dx 
£ 


Substituting the rate u 


6 u 

fit '3 £V 




into (A. 3) and requiring that the energy flow ft is converted into the surface 
energy 


we obtain 


SE = 4£y = 2£G , G c = specific fracture energy 

as defined by Irwin 


a 

2Y J [( 1^ + • fldx = G c 

£ 


(A. 5) 


(A. 6) 


* Note that in the subcritical range load s' is treated as a function of crack 

length. The function 6"= & (t) is a priori unknown and will be subject to determination. 
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This can be further reduced to (Wnuk [6]): 

a 

£ 

which indeed is the equation governing the slow growth in the subcritical 
range. For a rate- dependent solid we shall have to restrict the validity of the 
above governing equation to small scale yielding, since only then the viscoelastic 
displacement u(x,o,£) and its increment fiu can be expressed as a product of the 
elastic solution to a given boundary value problem u° = u°(x,a,£) and the 'creep 
function" ip = ip ( t) , as follows (Wnuk [6]): 


u(x,o,H) = u°(x,o , £) ■ ip(6t = A /£) , A = material constant 

5u(x,a,£) = fiu°(x,a,£) • ip( 6 1) + u° fiip - 5u°(x,a,£) • ip(6t) 


(A. 8) 


Note that the argument of the creep function ip has been replaced by the time 
interval fit equal roughly the time used by the crack tip to traverse its own 
plastic zone. The essential assumptions made here are that the time interval fit « A/ft 
is sufficiently small and that the function ip does not vary rapidly within the in- 
terval fit. This of course implies small scale yielding (A -*■ 0) range, while the 
rate of loading should be restricted to the "slow" one. These conditions are 
satisfied when the "inherent opening distance" A is much smaller than the character- 
istic length £* = ttK 2 /8Y 2 , and when an increment of load Aa associated with time 
equal to the representative relaxation time is negligible. With these assumptions 
equation (A. 7) simplifies considerably, and it reads 

a 

2Y §irJ u ° (x >°>* )dx + = why (A - 9) 

Z 


Here only the knowledge of elastic solution is implied, and thus both terms in eq. 
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(A. 9) can be identified with the linear fracture mechanics entities, namely 


do" 5 

2Y | u°(x,o,£)dx = M = slew growth operator 


(A. 10) 


£ 


and 

2Yu°(tip) = G = Irwin's energy release rate (A. 11) 

The first term M has been recently introduced into the theory of subcritical propa- 
gation (reference 6) , and may be thought of as a measure of inelastic behavior of 
a solid. The other term is a well-known Irwin’s energy release rate G or Rice’s 
path independent integral J (or the derivative of the strain energy G - -^(dU/dZ) J. 

It may be readily verified that the Grif fith-Irwin instability point is included 
in the equation (A. 9), but it results from it only 

i) if the solid is perfectly elastic, then M = 0 and = 1. 

ii) when the final point in the succession of meta-stable states 

d0 

(i.e. for a slcwly growing crack) is attained, as then ~ = 0, and again one re- 
covers Irwin's equation G = G c , provided that the matrix is rate-independent. How- 
ever, if the rate-dependence is there, the simplest form of eq. of motion (A. 9) 
follows for a crack growing under a sustained constant load (do =0), and it is 

G(o , £) = G c /iK A/£) (A. 12) 

We shall briefly illustrate how the above equations, in particular (A. 9) and 
(A. 12) work. To do so we shall expand the function ip (<5 1 ) in the Taylor series 
around the point 6t = 0 (implying smallness of the interval 6t = A /£ , which is 
satisfied for £ > 0, that is for a crack which already moves. Applied load, there?- 
fore, must be above the threshold value). We have 


6 1=0 


*<«t - o) - i + [■fSM] 


• St + . . . 


(A. 13) 
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or ? denoting the derivative [df(t)/dt] = B and recalling that 

1 t=0 


we get 



A_ da 
d di 


*(«t) = 1 + B ^ 


(A. 14) 


(A. 15) 


for an arbitrary linear viscoelastic solid. Finally computing M (see reference 6) 
we reduce the governing equation (A. 9) to the form 


G(o,£){ 


l + ^W . 3g [1 + B 442- ] = 


12nY‘ 


3a 


a dH dH 


l+B 

a d£ 


(A. 16) 


in where n equals unity for plane stress and 1-v 2 for plane strain. 

The above equation can be linearized if we restrict our attention to the near- 

2 3 

critical states only and assume that the powers (~~) and (— — ) are negligible 

Q Q-X/ 

vs. Then (A. 16) reduces to 


1 + t 


E + ^ f = G c /G<a,£) 


(A. 17) 


Consider a plane crack contained in a large plate subject to uniform tension, then 


6(a,£) = na 2 (tt£) /E (a) (A. 18) 

Second configuration of interest is that of a plane crack opened by a pair of 
point forces applied directly to the crack faces, then 

6(o, £) = (b) (A. 19) 


For these two configurations eq. (A. 17) takes the form 


(— cB + C/B) — = 2 

' dc cB 2 


(a) 


f 2 1 + c/e) dB . gnji- ,B 2 

[3TT 2 z dC' B 2 


(b) 


(A. 20) 
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where the dimensionless load 3 and the dimensionless crack length £ are defined 
as follows 


6 = 


TTQ 

2Y 

ttP 

2Y£* 


for 

for 


(a) 

(b) 



both cases 


Two limit cases which can be treated somewhat easier than the general case are 
worth mentioning here, and these are 

1) [3 = 0, constant load fracture (static fatigue), for which the equation of 
motion can be integrated in a closed form. 

2) In the other limiting case (3 00 , fast loading or time- insensitive matrix) the 

governing equation simplifies again, since now C/S = 0. Fig$J a and b show the 
family of integral curves resulting from eq. (A. 20a). They include two extremes 

of possible behavior: 

a) fast loading, cf. the steepest curve; 

b) sustained load (zero rate); cf. the horizontal line. 

This clearly demonstrates that, all other factors being fixed, the rate of loading 
and the initial crack length determine the critical stress at which failure will 
occur. In a similar way the delay time, which elapses between load application 
and the final instability, is affected by the rate-sensitivity of the solid. 


TABLE II 


Recorded numbers of cycles to failure at various stress 
ranges A and sensitivities £. The first row gives 
the number of cycles between J = 1 and 7,-10, the 
second one refers to 10 < g < 100, while the 
third row gives the total numbers of 
cycles to failure. 


|ru/» 

Afl 

^ max 

RATE SENSITIVITY 


C- 1.5 

1.0 

0.5 

0 



5587 

7876 

13248 

47082 

.03 

.13 

2084 

2415 

2884 

3978 



7673 

10291 

16132 

51060 



3771 

5328 

8944 

30999 

.04 

.14 

1365 

1585 

1902 

2482 


5146 

6913 

10846 

33481 



2713 

3846 

6457 

21491 

.05 

.15 

937 

1088 

1311 

1713 



3650 

4934 

7768 

23204 



2019 

2890 

4861 

17966 

.06 

.16 

676 

782 

943 

1245 



2695 

3672 

5804 

19211 



1563 

2250 

3728 

11978 

.07 

.17 

498 

568 

694 

920 



2061 

2818 

4422 

12898 



1237 

1775 

2994 

9390 

.08 

.18 

363 

424 

520 

689 


1600 

2199 

3514 

10079 



998 

1426 

2420 

7480 

.09 

.19 

280 

323 

408 

525 



1278 

1749 

2828 

8005 



822 

1168 

1992 

6062 


.20 

218 

255 

309 

409 

. 1 

1040 

1423 

2301 

6471 




TABLE III 


Recorded numbers of cycles to failure at various mean stresses 
0 me<an and sensitivities C . The first row gives the number 
of cycles between^ = 1 and t, - 10, the second one 
refers to 10< J < 100, while the third row 
gives the total numbers of cycles 
to failure. 


A 3 = .01 


f 



Stress 

level 


c = 0.0 

C = 1 


io 

II 

ro 

.12 

4 (3 

6 .13 


107886 

8742 

19642 

5541 


10187 

4192 

( 

TOTAL 

116628 

25183 


14379 

.13 

< ft 

£ .14 


85306 

6989 

16539 

4454 


8400 

3344 

" r 

TOTAL 

92295 

20993 


11744 

.14 

£ (3 

^ .15 


70533 

5553 

13819 

3567 


7112 

2625 

\ 


TOTAL 

76086 

17386 


9737 

.15 

] 

< ft 

^ . 16 


58475 

4350 

11460 

2794 


6362 

2072 

( 

TOTAL 

62825 

14254 


8434 

1 

. 16 

£ f 

* ^ .17 


49769 

3406 

9862 

2176 


5349 

1655 

V 

TOTAL 

53175 

12038 


7004 

.17 

£ £ 

$ £ .18 


42351 

2794 

8579 

1899 


4732 

1388 

l 

TOTAL 

45145 

10478 


6120 

.18 

- P 

£ .19 


36367 

2263 

7444 

1492 


4164 

1157 

l 

TOTAL 

38630 

8936 


5321 

.19 

<. i C 

5 £ .20 


32669 

1875 

6779 

1278 


3888 

983 

v r 

TOTAL 

34544 

8057 


4871 




TABLE IV 


Testing the cumulative damage law on the EAI380 analogue computer. 
Note that the deviation A increases for greater 
sensitivities and discrepancies in 
the applied stress level. 


^min'P'lfVax 
.19 - .20 

.12 - .13 


32669 = 946 
34544 


8742 


116628 


- .075 


total = 1.021 


C = 1 
6ZI2 = .841 


8057 
4454 _ 


20993 


- .220 


total = 1.061 


C = 2 

3888 — # 79g 
4871 


3344 _ 


11744 


- .292 


total = 1 .090 


.18 - .19 
.13 - .14 


36367 _ 
38630 


7444 „ 

8936 


= .833 


4164 _ 


6989 _ 


92295 


.076 


4454 _ 


20993 


.212 


5321 
3344 _ 


- .783 


11744 


= .285 


total - 1.017 


total - 1.045 


total - 1 .068 


.17 - .18 
.14 - .15 


12351 = .938 

45145 

555 3 = >Q73 
76086 


8579 _ 


10478 
3567 _ 


17386 


.819 


4/42 _ <773 
: 6120 


CM 


.205 


i 2625 „ .j 7 n 

+ 

9737 


total — 1.011 


total = 1 .024 


total =1.043 


16 - -17 
,15 - .16 


REVERSED 

ORDER 


19Z69 = .936 
53175 


98 6 4 - .819 
12038 


5212 _ _ 764 

7004 


4350 


62825 


= .069 


2794 


14254 


= -196 


2072 


8434 


= .246 


total = 1 .005 


total = 1-015 


total = 1.010 


Ni y N i 

— — i r? — 


Ni 


if 


Na 


if 


8 X 


■o 


Nj f y_ Nj 


'if 


Ni 


if 


o » 

+5 fH 

« c 
■ r( K 

>'. 

Q) - 


N if 


7 N i 

) 4 — • ^ — 
i N if 


£ 

8*1 

o w 
io® 

•H C 
>: 

TJ 


,12 - .13 


.19 - .20- 


107886 _ q „ 
116628 


19642 _ 


1875 

34544 


= .054 


25183 
1278 _ 


= .780 


10187 _ 


8057 


= .159 


14379 

983 


= .708 


= .202 


4871 


total = 0.979 


total = 0.939 


total = 0.910 




13 - .14 


18 - .19 


85306 


92295 


o= .924 


16539 


R 263 . _ 059 
38630 


20993 

1492 


8936 


- .788 


= .167 


8400 


11744 

1157 


5321 


- .715 


- .217 


total = 0.983 


total = 0.955 


total =0.932 


14 - .15 


17 - .18 


70533 _ 


76086 


= .927 


13819 = i?95 


7112 


2Z2^ = .062 

45145 


total = 0.989 


17386 

Mii-= .181 

10478 


9737 

1388 


6120 


= .730 


= -227 


total = 0.976 


total = 0.957 


,15 

,16 


.16 

.17 


58475 


53175 


= .931 


= .064 


11460 _ 

14254 “ *" 04 


6362 _ 


2176 


12038 


= .181 


8434 

1655 

7004 


.754 

.236 


total = 0.995 


total = 0.985 


total =0.990 
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Fig. lb. Subcritical growth induced by a monotonic load at three different 
rate sensitivities; simulated by the analogue computer EAI 380. 
Note that the increased internal friction (C)enhances the amount 

of slew growth. Arrows indicate points of terminal instability. 
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FIG. 3 EFFECT OF THE INITIAL FLAW SIZE CW AND THE RATE SENSITIVITY (C) 
ON THE AMOUNT OF SLOW GROWTH PRECEDING CATASTROPHIC FAILURE 

( ACCORDING TO DATA GENERATED ON EAI 380 ANALOGUE COMPUTOR) 




Fig. 4a. Analog computer diagram for integrating 


m s (t)x. 

/ - m*(t) x. 


uation X — C 
































Fig. 5a, Effect of the rate sensitivity C_ = C/< on the fatigue crack growth. 

Increased^ is equivalent to enhanced sensitivity or/and lower frequency. 

Here: C = 0, initial crack length j>®= 1 

stress range 0.1 ^ ^ 0.3 

number of cycles not shown = 169 

total number of cycles = 426 






(b) C = 0.5, initial crack length Jj* = 
stress range 0.1 ^ {b $0.3 

number of cycles not shown = 160 
total number of cycles = 300 


(c) = 2.5, initial crack length 1 

stress range 0.1 p $0.3 
number of cycles not shown = 92 

total number of cycles = 132 


Fig. 5b & c. Effect of rate sensitivity = C/<[(i> on the fatigue crack , grow th . 

Increased C is equivalent to enhanced sensitivity or/and lower frequency. 
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Fig. 6 Final stage of fatigue life simulated by the EAI 380 analogue computer. 
Note the enhanced rate of growth at larger value of sensitivity C. 





Fig. 7 Effect of stress range p on the fatigue life (see Table II for the 
measured numbers of cycles to failure). 






Fig. 7d. Photograph from the oscilloscope screen of runs 7a, 7b 
and 7c combined. 
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Fig. 8 Effect of the mean stress on the fatigue life (see Table III 
for the measured numbers of cycles to failure) . 
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Fig. 11 Subcritical growth of a crack under (a) uniform tensile field, (b) point 
forces applied directly to crack surface; (c) shows loads imposed, 
reduced load, crack length and time are; m = 6"(t) f6 Ck . r , x = Jt/ 1 ,-0=t i^(0) 





